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Abstract 

Effective theories of a scalar f invariant under the internal galileon symmetry 
f —>• <f> + b^x^ have been extensively studied due to their special theoretical and phe¬ 
nomenological properties. In this paper, we introduce the notion of weakly broken 
galileon invariance, which characterizes the unique class of couplings of such theories 
to gravity that maximally retain their defining symmetry. The curved-space remnant 
of the galileon’s quantum properties allows to construct (quasi) de Sitter backgrounds 
largely insensitive to loop corrections. We exploit this fact to build novel cosmolog¬ 
ical models with interesting phenomenology, relevant for both inflation and late-time 
acceleration of the universe. 
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1 Introduction and summary 


Over the years, the effective field theory (EFT) approach has proven itself invaluable for or¬ 
ganizing physics at different length scales. With the sole knowledge of the relevant physical 
degrees of freedom and symmetries - exact or approximate - that govern their dynamics, 
EFTs can be used to infer model-independent physical predictions. The reason for such uni¬ 
versality lies in a fundamental property of the great majority of physical systems, according 
to which the details of physics at short distances do not leave a qualitative impact on their 
large-distance characteristics, enabling one to analyze theories with various UV structure in 
a single framework. 

In recent years, the EFT approach has offered new insights in various cosmological 
contexts, including - but not limited to - inflation in ei, dark energy mm and the large- 
scale structure [6] . For example, the effective theory of inflation is based on the observation 
that the dynamics of the most general theory of ‘single-clock’ inflation can be universally 
captured by an EFT nonlinearly realizing time diffeomorphisms (cliffs (Q t — > £ + £°(i, x) , with 
spatial cliffs x l — > + x) realized linearly. The energy scales of interest are those around 

the inflationary Hubble rate, E ~ H, i.e. the frequency at which all correlation functions 
are measured in the CMB. The spectrum of perturbations consists of the two polarizations 
of the graviton plus the Goldstone boson of time translation symmetry breaking. That the 
latter mode has to be present in the spectrum is a direct consequence of the symmetry 
breaking pattern and has little to do with the exact UV details of the microscopic theory 
of inflation. This Goldstone boson is what we usually refer to as the adiabatic mode. The 
CMB indicates that precisely this mode is predominantly responsible for the generation of 
large-scale structures in the universe. Indeed, any UV theory of inflation that does not lead 
to extra degrees of freedom around the Hubble energies is equivalent to single-clock inflation 
and is thus captured by the EFT of Refs. DEI, 

Given that symmetries define the effective theories, it is of interest to explore the 
possible symmetries of systems consisting of one or more scalar fields coupled to gravity - 
a typical setup for cosmological model-building. For example, in the context of inflation, 
effective theories with (approximate) shift symmetry play a key role. In this paper, we 
study yet another possible - and, as we argue below, necessarily approximate - symmetry 
of cosmological scalar fields: the invariance under internal galileon transformations 


->■ 0 + b 


( 1 . 1 ) 


Theories invariant under (1.1) have appeared in various contexts before. To start with, 


Eq. (1.1) is a symmetry (up to a total derivative) of the simplest possible quantum field 
theory: that of a free scalar field. The simplest interacting generalization, nontrivially 


invariant under (1.1), i.e. containing less than two derivatives per field in the Lagrangian, 


Hn what follows, we will at times abuse the nomenclature by referring to this redundancy as ‘time 
translations’. 






has appeared in Ref. [7j in the context of the DGP model [Hj. The most general scalar 
theory with the same property - the galileon - has been proposed in [9] and subsequently 
found [10] to describe the scalar polarization of the ghost-free dRGT massive graviton HU- 
Galileon theories are interesting in many ways. Their defining property is that, despite being 
higher-derivative theories, they have no more than second-order equations of motion, thus 
describing a single propagating scalar mode. Moreover, as a direct consequence of invariance 


under (1.1), the coefficients of the leading galileon interactions are not renormalized, at least 


perturbatively [7], making any tuning of them stable under quantum corrections. 

Despite these interesting properties, the invariance under galileon transformations can¬ 
not be exact in nature because the couplings of the scalar to gravity, at least, break it 
explicitly. We will therefore be interested in characterizing the theories that preserve as 
much as possible the attractive quantum non-renormalization properties of the exactly in¬ 
variant case. This naturally leads to the notion of weakly broken invariance under (1.1), 


which we will define in a precise way in what follows. The simplest sufficiently non-trivial 
theory with weakly broken galileon invariance (and exact shift symmetry) is of the following 
form 


c = -)(a« 2 + 4w) 2 n^ + Aw) 4 


At 




( 1 . 2 ) 


While the first two operators in (1.2) are exactly invariant under ( 1.1[ ) (up to the boundary 
terms), the quartic operator is a small breaking as far as A 2 A 3 . In general we would 
expect other symmetry breaking operators of the form (<9</>) 2n to be generated by quantum 
corrections at the scale A 2 . However, in this case, a stronger result, the remnant of the 
non-renormalization properties of the invariant action, holds: all the symmetry breaking 
operators are generated at a scale that is parametrically higher than A 2 . This means in 
particular that the operator (<90) 4 gets only small corrections through loop effects. In the 
presence of gravity, as we will show in the next section, the weak breaking gives rise to even 
more non-trivial consequences. 


Importantly, the theory (1.2) admits a large generalization that retains both its quan¬ 


tum properties and the second-order equations of motion. Perhaps not surprisingly, we will 
find that the proper generalization fits into the class of the most general Lorentz-invariant 
theories with second-order held equations, known as Horndeski theories CGI- Since, however, 


a generic Horndeski theory does not have much to do with the invariance under (1.1), we 


prefer to refer to the subclass that we study here as ‘theories with weakly broken galileon 
(WBG) invariance’. The purpose of this work is to introduce these theories and to set the 
stage for their detailed phenomenological studies. 

The study of galileon-invariant theories in the inflationary context have been initiated 
in Ref. M- The model, referred to as ‘galileon inflation’, is based on the particular curved- 
space extension of the theory that keeps the property of second-order equations of motion 
intact - the so-called ‘covariant galileon’ of Ref. [13]. It has been pointed out in [13], that 
such theories enjoy more freedom compared to ghost/DBI inflation-like models as far as 
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phenomenology is concerned. In particular, they can lead to the possibility of lifting the 
‘large speed of sound/small non-gaussianities’ correspondence, characteristic of the latter 
models. Our results extend the findings of Ref. ra in several directions. While it has been 
assumed in [13] that the covariant galileon is the unique model capable of extending the 
phenomenological virtues of shift-symmetric theories in a radiatively stable way, we find that 
there is in fact a wider class of models that can achieve this. Unlike the covariant galileon, 
the theories we will be interested in do not generically reduce to the standard galileon 
once gravity is turned off. Nevertheless, neither non-renormalization nor the second-order 
equations of motion of the galileon need to be given up, leading to the possibility of strong 
- and quantum-mechanically robust - phenomenological differences from slow-roll inflation. 
This generalizes the well-studied case of DBI models in an interesting way. 

The paper is organized as follows. We start in Secs. [ 2 ] and [3] by showing that a (shift- 
symmetric) subset of Horndeski theories can be derived solely based on the requirement 
of WBG invariance at the quantum level, a concept that we define along the way. In 
Sec. [4] we study general cosmological implications of approximate symmetry under galileon 
transformations and apply our findings to inflationary, as well as late-time cosmologies in 
Secs. [5] and [6j respectively. Finally, in Sec. [7] we conclude. Various technical calculations 
that would overwhelm the main body of the text are collected in the two appendices. 


2 Galileon invariance (with and without gravity) 


In this section we study the fate of galileon symmetry upon inclusion of gravity, starting out 


by reviewing flat-space theories exactly invariant under (1.1). If internal galileon invariance 


is realized in nature, however, it has to be approximate: the coupling to gravity unavoidably 
breaks this symmetry and, at the same time, defines the sense in which the breaking can 
be considered weak. Indeed, as we will see below, out of the infinite number of possible 
ways one can couple the galileon to gravity, there is a unique set of non-minimal couplings 
that can be qualified as being ‘more invariant’ under galileon transformations than all the 
rest. The resulting theory is the covariant galileon of Ref. [H]. On the other hand, since 
the galileon symmetry has to be weakly broken by the couplings to gravity, one should in 
principle include in the effective theory symmetry-breaking operators of e.g. the form ( d(j)) 2n . 
The latter will be generated by quantum loops with Wilson coefficients bounded from above, 
so as to be consistent with the approximate symmetries at hand. As we will see, requiring 


invariance under (1.1) to be only weakly broken will naturally lead to a particular sub-class 


of Horndeski theories, which is however much more general than just the covariant galileon. 
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2.1 Flat space galileons 


Consider a trivial, free theory of a scalar </>. As emphasized above, in addition to more 
familiar symmetries (such as the ones under constant shifts or conformal transformations), 


this theory possesses an extra invariance under internal galileon transformations (1.1). The 
latter leaves the action invariant only up to a boundary term. As shown in [2j, apart 
from the free theory - and even a more trivial tadpole term, there are exactly three more 
interaction terms in four spacetime dimensions that share invariance under internal galileon 
transformations in a non-trivial way. The corresponding theory can be written as 


£ = w) 2 +£ AM. 

1=3 A 3 


with the three interaction terms given by 


£3 = W) 2 [$] . 

A = W) 2 ([4f-[4> 2 ]) , 

£ 5 = (a« 2 ([t] 3 -3(®][* 2 ]+2[« 3 ]) 


( 2 . 1 ) 


( 2 . 2 ) 

(2.3) 

(2.4) 


where we denote [<&] = D(p 1 [4> 2 ] = etc. In addition to being invariant under 

(1.1), galileon theories share another special property: the associated scalar equations of 


motion are second order, both in time and in space, despite the higher-derivative interactions 
in Eqs. (2.2)—(2.4). This guarantees that there are no extra propagating degrees of freedom 


hidden in </>. Moreover, the particular structure of the Lagrangia.n (2.1) results in a non¬ 


renormalization theorem that allows to radiatively generate only terms trivially invariant 


under (1.1), i.e. terms with at least two derivatives acting on (ft. In particular, the coefficients 


C/ of galileon interactions are not renormalized by quantum loops in the presence of exact 
galileon invariance [7j. In theories of modified gravity with matter, galileons couple to the 
rest of the degrees of freedom via symmetry-breaking Planck-suppressed operator^, usually 
making the renormalization of C/ parametrically weak. 


2.2 Coupling to gravity 

It is generally impossible to couple the galileon to gravitjj^J while keeping invariance under 
galileon transformations, or their curved-space generalization, intact. Thus, galileon sym¬ 
metry is expected to be broken even in the purely scalar sector - i.e. if one sets the metric 
perturbation to zero ‘by hand’ - due to loop-generated operators not invariant under fll.l ) 
(these operators of course have to be suppressed by at least one power of the Planck mass). 

2 In massive gravity, where <j> describes the helicity-0 polarization of the graviton, this simply follows from 
the equivalence principle. 

3 Unless gravity is massive, see [15; . 


4 













d m 4 ), m> 2 


d n h 

'N/X/N/X/X/X/V 


90 


Figure 1: For every given vertex, a solid straight line corresponds to a possible external leg resulting 
in a scalar with at most one derivative acting on it. A dashed line corresponds to an external scalar 
with at least two derivatives, while a wiggly line denotes an external graviton. 


For the sake of concreteness, let us concentrate on the operators of the form (<90) 2n . After all, 
it is the absence, or smallness, of these operators that makes the galileon special compared 
to a generic shift-symmetric theory^] These operators can be generated from loops with the 
appropriate number of insertions of symmetry-breaking vertices, that include at least one 
graviton line. Interpreting A 3 - the smallest scale by which interactions are suppressed in 
the effective theory^ - as the genuine cutoff, any loop-generated operator can be written as 


(90) 2n 

A 4(n—1) 

IV k,n 


A*,„= [MfiAf-*- 4 ] 4 -- 


(2.5) 


where k denotes some positive integer and we have omitted factors of 167 T 2 for simplicity. For 
fixed k and n sufficiently large, the scale suppressing a given single-derivative operator, Afc jn , 
can in principle be arbitrarily close to A 3 . If this were true, the resulting theory would by no 
means be considered as a theory with WBG invariance, since symmetry-breaking operators 
would be order-one in the units of the cutoff. 


We wish to show here that: i) this in fact does not happen for the galileon, even 
if minimally coupled to gravity; ii) for a very particular, non-minimal coupling, the scale 
suppressing the symmetry-breaking operators can be made parametrically higher than the 
analogous scale characterizing all the other ways of introducing gravity into the system. 
For the latter theories, there is a well-defined separation between the scale suppressing 
the invariant galileon interactions and the quantum-mechanically generated single-derivative 
operators: while the former are suppressed by powers of A 3 , the scale suppressing the latter 
is parametrically higher, asymptoting to 

A 2 = (MpAj ) 1 ^ 4 5 (2.6) 


for large n in (2.5). Perhaps not surprisingly, we will fold that the special type of couplings 
to gravity required to suppress symmetry-breaking operators is of the Horndeski class - 
the unique curved-space extension of the (generalized) galileons that leads to second-order 
equations of motion both for the scalar and the metric. 


Let us now prove all of the above-described statements, starting with a few definitions. 
Consider an arbitrary vertex of the form (<90) fc (<9 m 0) n 9 / /i p , with m > 2. In terms of Feynman 

4 Moreover, the size of these operators is an avatar for technical naturalness of non-trivial de Sitter vacua 
of the classical theory, as we will show in the next section. 

5 We always assume in this paper that A 3 <C Mp\. 
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diagrams, we will adopt the convention according to which a solid straight line corresponds to 
an external leg coming out of this vertex, which results in a scalar with one derivative acting 
on it. A dashed line corresponds to an external scalar with at least two derivatives, while 
a wiggly line denotes an exernal graviton, see Fig. [l] For example, a fc = 2, n — 1, p — 1 
vertex would correspond to the first diagram on the upper line of Fig. [2j In certain cases, 
the number of solid lines can be less than the number of factors of <90 in a Lagrangian 
interaction term. For example, the special structure of the pure galileon interactions makes 
them equivalent to vertices with only dashed lines, even though on average these terms 
contain less than two derivatives per field. This is because a galileon vertex can only lead 
to external states with at least two spacetime derivatives - the fact that lies at the heart of 
the non-renormalization theorem associated with these theories [7J. We now wish to show 
that there exist more redundant vertices of this kind in the suitable non-minimal extension 
of the theory to curved space. 


Let us have a look at all possible vertices with a single graviton line in the minimally 
coupled galileon theory, obtained from (2.1) by replacing all derivatives with covariant ones 
(<9 M —> V /( ). A straightforward inspection yields that there are six of such vertices, shown in 
Fig.! Vertices with three solid lines can in principle arise from picking up a factor of dh 
from covariant derivative^ i.e. V 2 0 ~ d 2 (f) + dhdcp, where we denote by V 2 any product of 
two covariant derivatives applied on 0. This means that, in the minimally coupled galileon 
theory, the smallest scale by which the operators of the form (2.5) with n = 3a, k = 
2a are suppressed, T r\j (MpiAl) 1 / 6 . The latter suppression scale is still too small to be 
consistent with the generic definition of WBG invariance, which we will introduce and discuss 
extensively in what follows. 


We will now show that it is possible to modify the theory by non-minimal couplings 
to gravity that result in the elimination of vertices with three solid lines, leaving just a 
factor of (<90) 2 per graviton line. Insertion of these vertices into a generic loop diagram 
can only lead to symmetry-breaking operators of the form (2.5) suppressed by at least the 
scale A 2 = (MpiAg) 1 / 4 ^ parametrically higher than (MpiA|) 1//6 . Therefore, there is a well- 
defined sense in which the resulting theory is ‘more invariant’ under the internal galileon 
transformations than a generic curved-space extension of the galileon. This defines what we 
mean by ‘theories with WBG invariance’ throughout the present work and we will sharpen 
this definition in the next section, where non-trivial vacua of the resultant theories are 
considered. 


For the cubic covariant galileon, it is easy to show that the first vertex from the lower 
line of Fig. [2] is absent: the only way to generate it is by picking up a metric perturbation 
from the covariant laplacian. This yields a term of the form (<90) 3 <9h. Putting the derivative 

6 Moving the derivative on h to the rest of the fields in the vertex generically does not reduce the number 
of solid lines for the quartic and quintic galileons, since this extra derivative can act on the factor d 2 (f>, 
representing the dashed line. 

'This can be seen by inserting enough number of vertices with one graviton and three solid lines into a 
generic 1PI loop diagram. 
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Figure 2: Single-graviton vertices, which can contribute two or three external scalars with one 
derivative (( dct>) 2 or ( dq i>) 3 ) to a 1PI vertex. The vertices from the second line exactly cancel for the 
unique (non-minimal) curved-space extension of the galileon that retains second-order equations of 
motion. 


from h on the rest of the holds via partial integration however makes it evident that the 
corresponding vertex can only have two solid lines, but not three. The case of quartic 
and quintic galileons is more non-trivial, but straightforward; we show in Appendix |A| that 
vertices with three solid lines and one graviton, as well as five solid lines and two gravitons 
can be removed by suitably adding non-minimal couplings to gravity. The resultant theory 
is the covariant galileon of Ref. [TJ] — characterized, as a bonus, by second-order equations 
of motion both for the scalar and the metric. 


In summary, the ‘most symmetric’ generalization of the galileon coupled to gravity 
consists of the following operators 


1 _ 1 

AT 3 A| 

1 _ 1 

Ap 4 V ~ 9 Af 

Xf 5 ^ Ajj 


/’min 
,3 ? 


(d<l>)‘ if - 4 £f ” , 

W) 1 + ? Cf 


(2.7) 

( 2 . 8 ) 

(2.9) 


where Cf ia denote the galileons (2.2)—(2.4) minimally coupled to gravity. The structure of 
the full effective theory is such that every pair of external scalars with no more than one 
derivative on each unavoidably comes with a suppression of at least one power of the Planck 
scale. With this in mind, quantum-mechanically generated operators of the form (90) 2ri can 
be estimated simply on dimensional grounds to be at most of the following magnitude 


(90) 2n 

MpjAp 4 ' 


( 2 . 10 ) 
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The scales (Mp L .Aj n 4 ) 1 /( 4n 4 ) suppressing such operators approach A 2 only for asymptotically 
large n, otherwise being parametrically larger. 


In the following sections, we will argue that the statement of approximate galileon 
invariance can be viewed as a statement about non-trivial classical vacua, generically present 
in the shift-invariant theories at hand. If these vacua are to be insensitive to UV physics, 
the symmetry-breaking operators in the effective theory can be at most of order (2.10) in 
magnitude. In such cases we will say that the galileon invariance is broken only weakly by 
couplings to gravity. 


3 Theories with WBG invariance 


Perhaps, the most important lesson that one can draw from the discussion of the previous 
section is that the covariant galileon is in fact not the most general theory enjoying the above- 
described quantum properties. Consider the effective theory with the leading operators given 
by the following subclass of the Horndeski terms 


£r G = G 2 (X) , 

•cr° = g 3 (x)m , 


A! 


/’WBG 

^4 


A 8 A 4 

= xi G,(X)R + 2jl G lx (X) ([4f - [« 2 ]) , 


A§ 


/’WBG 

J ~'5 


A 8 A 4 

X§ G s (X)G^ - ^ G«(I) ([*] 3 - 3[4>] [« 2 ] + 2[« 3 ]) 


(3.1) 

(3.2) 

(3.3) 

(3.4) 


where we now extend the notation employed in Eqs. (2.2)—(2.4) by replacing the partial 


derivatives with covariant ones, i.e. [<f>] = g^X u (j), [<f> 2 ] = V M 0, etc. Moreover, 

Gj are arbitrary dimensionless functions of the dimensionless variable 


X = g^d^dvcj) , 

and the subscript ‘A"’ means differentiation with respect to this variab 1 ^ 
assume that the functions Gj are defined through the Taylor expansion 


Gr = c 


( 0 ) 


c^A + c y ;’X z + • ■ • = cr + C 


„( 2 ) 


(0) 


M) 


m 2 

A 4 


+ c 


( 2 ) 


A 8 


+ .. 


(3.5) 

We will then 

. (3.6) 


8 The particular form of the interactions (3.1 )~( ph4| ) (the relative coefficients between two operators in 

/’WBG or £WBG^ f or example) may appear tuned, and this is sometimes presented as an unfortunate feature 

in the literature. We stress that this ‘tuning’, motivated by restoring unitarity in the theory, is in fact natural 

and stable under quantum corrections. This is very similar to why we work with a gauge-invariant kinetic 
term —Tr F^ 11 F^ u in theories with massive spin-1 particles, despite of there being no gauge invariance in the 
presence of the mass term. 










where c\ n ' 1 are dimensionless, order-one coefficients. Note that by setting all the coefficients 


cj"' ) to zero except c£\ c^\ cf ' 1 and Cg 2 ' we recover the covariant galileon of Eqs. (2.T)—(2.9). 


As in the example of the introductory section, the theory at hand is characterized by 
two scalc^ A 2 and A 3 <C A 2 . From the EFT standpoint that we are adopting here, all 
Wilson coefficients g\"' > are measured in the units of the cutoff A 3 . The theory (3.1)—(3.4) 
is the one for which a well-defined class of (canonically normalized) operators - those that 
respect the galileon symmetry 0 —)• 0 + - have order-one coefficients, 


(n) 

91 


Oil) 


(3.7) 


while all the others are parametrically suppressed with respect to the coefficients by 
positive powers of the ratio A 3 /Mpi. As a consequence, one can define a scaling limit 


Mpi —y oo, A 2 —y oo, 


A 3 = hnite 


(3.8) 


in which the system recovers the exact internal galileon symmetry. 

A comment about the terms and c^ is in order herej^j These terms are absent in 

the covariant galileon, which starts with G 4 , G 5 ~ (90) 4 . Thus, one can wonder whether they 
modify our conclusions regarding quantum corrections. By expanding £4 and £5 at linear 
order in metric perturbations, it is straightforward to verify that and Cg 1 " 1 do lead to new 
operators at the scale A 3 . However, these operators are just a scalar-tensor generalization of 
the galileon familiar from the decoupling limit of dRGT massive gravity [TO] . 


£ 


WBG 

4 

£WBG 


h e 

1 c 11.1;^ /, 




rs - / hp,v£papX^isPcrddoidpcftdpda-CfodxdfiCj) 4 “ • • • • 


(3.9) 


The latter interactions obey the same non-renormalization theorem as the galileon [T 6 ] , 
forbidding asymptotic states with less than two derivatives acting on them. In particular, 
this means that the counting of Mpj 1 factors in loop diagrams of the previous section goes 
through unaltered: each factor of (<90 ) 2 arising from a generic 1PI vertex still comes with a 
suppression of at least a factor of MG 1 , and our analysis of quantum corrections performed 


for the covariant galileon also applies for the generalized theory (3.1)—(3. 10 Therefore, for 
simplicity we will disregard these terms and assume that G 4 and G 5 start at least quadratic 
in A". Indeed, our discussion of the previous chapter guarantees that, once tuned to zero, G^x 
and G 5 x do not receive any appreciable quantum corrections as a result of the approximate 
galileon symmetry. 


9 The Planck scale arises here from £^ BG in Eq. (3.3) as A/p, = 2^ h- 2 / 


10 The operator c ^ is just the cosmological constant, which for simplicity we tune to zero, c ^ = 0. More¬ 
over, the operators c!j 1 and set the normalization of the scalar and graviton kinetic terms respectively, 


while Cg°^ and c^ UJ give inconsequential total derivatives and can be disregarded. 


JO) 


11 That the terms 
has been noticed in 


l3 1 9h can be rewritten as a certain limit of non-minimally coupled scalar-tensor theory 
117] and the cosmology of that theory has been studied in [18] . 
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Notice that the theories that we propose do not reduce to the galileon for general 
functions Gj, even if one switches off gravity. Rather, they describe a certain generalization 
thereof, which still has no more than second-order equations of motion. The symmetry 


(1.1) is thus broken even when one sets = 0 (note that this is different from taking the 


decoupling limit (3.8)). Nevertheless, the corresponding breaking of invariance under (1.1) 


is in a well-defined sense ‘not stronger’ than the one already present once gravity is turned 


on. In particular, it follows from the very construction of the operators in (2.T)—(2.9) that 


each pair of factors of <90 still comes with a suppression of at least a factor of Mpi in the 


full effective theory. Thus, the Lagrangians (3.1)—(3.4) define the most general scalar-tensor 


theory consistent with our definition of WBG invariance. 

The peculiar quantum properties characteristic of the theories at hand put loop cor¬ 
rections under complete control for a broad class of physical backgrounds. This can be 
established by noting that loop-generated operators always have at least one extra factor 
of Mpj 1 , compared to their ‘tree-level’ counterparts. Indeed, working in the units in which 


A 3 = 1, each factor of h.\ becomes equivalent to Mpi, and the non-minimal terms in (3.1) 


(3.4) can be schematically written as 


MpiR(<90) 5 


i + ^! + ... 


Mi 


pi 


(V 2 0) r 


(3.10) 


while the similar loop-generated terms can never have the Mpi-enhancement. Likewise, the 
symmetry-breaking scalar operators generated quantum-mechanically are bound to have the 
following schematic form 


(90) 2n 
M n 

7Wpi 



(3.11) 


and the analogous terms in the original action are all enhanced by at least a factor of Mpi 
compared to this expression. Restoring the factors of A 3 , we conclude that if the held 
expectation values are such that 


1*1 £ 1 , 


\Y\ £ 1 


(3.12) 


where X is defined in Eq. (3.5) and 


Y = 


At 


(3.13) 


the loop corrections are negligible: they can never compete with the operators from the 


classical action (3.1)—(3.4) and all predictions based on the latter can be trusted in the full 


quantum theory. Moreover, not only are the magnitudes of the various operators technically 
natural, but any possible tuning of these can only be spoiled by corrections of higher order 
in A 3 /Mpi. This generalizes the galileon non-renormalization theorem in the presence of 
gravity. 
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On the other hand, if the conditions (3.12) are not satisfied, the background becomes 
too strong to be trusted: the magnitude of loop corrections becomes pumped up beyond that 
of the leading operators and galileon invariance alone is useless in organizing the low-energy 
expansion. One could only trust the classical theory beyond that point if there is an extra 
structure, leading to a resummation of the series (3.6). An example would be, e.g., the 
DBI theories [HE HI] (as well as generalizations thereof [21, 22j), where this happens due to 
non-linearly realized higher-dimensional spacetime symmetries. However, here we want to 
remain as general as possible and do not assume any extra symmetry beyond the (weakly 
broken) galileon one. In the latter case, the requirement of quantum stability expressed by 
Eq. (3.12) places a strong upper bound on how large the predictions for various physical 
quantities can be. 


4 De Sitter vacua 

Perhaps, the best way to illustrate the physical consequences of the discussion of the previous 
section is to resort to a concrete example. To this end, we wish to consider the following 
theory 


S 




v(4) + J2 4 VBG + • • • 

1=2 


(4.1) 


where CY BG are the operators defined in Eqs. (3.1)—(3.4) of the previous chapter. We have 
extracted and canonically normalized the scalar and graviton kinetic terms, so that G 2 , G 4 
and G , 5 are assumed to start at least quadratic in X (see the discussion below Eq. (3.9)), 
while Gs can have a linear piece. I 11 general, one can allow for a small potential, 


v(4>) 



(4.2) 


where the parameters A, m^, etc. can be (naturally) arbitrarily small, since they are the only 
ones that break the scalar shift symmetry (which otherwise is exact even on curved space). 
In the next section we will discuss the implications of a weakly broken galileon symmetry 
for inflation. We will thus be interested in background solutions given by an approximate de 
Sitter space with the Hubble parameter H. This section is devoted to a preliminary study 
of such solutions. 


For the theory (4.1), the scalar equation of motion on a flat FRW background reads 


1 d 
a 3 dt 


[2a 3 ^F(X,Z)\ 


dV 

dcj) 


(4.3) 
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where the function F is given by the following expression 12 

Gax 


F{X, Z) = - + G 2X ~ 3ZG 3X + 6 Z 2 [^ + 2G 4XX 


+ Z d 


G 


5X 


X 


+ 2 G 5XX ) , (4.4) 


with the variable Z dehned so as to roughly coincide in the limit (j) -C H<p with the back¬ 


ground value of Y dehned in Eq. (3.13), 


Z = 


Hep 

aT 


(4.5) 


Moreover, the two Friedmann equations can be written in the following way 

G 


3 M^H 2 = V + KX 


1 G2 , 9 r 

2 - T + 2G2A ' 

2 Z 3 I r,-+ + 2 G 


6 ZG: 


3X 




X 


4 G 


4XX 


X 


5XX 


(4.6) 


M^H = 


A \XF + - 4 ZG 4X - 8 ZXG 4XX - 3 Z 2 G 5X - 2 Z 2 XG 5XX ) 


1 + 2G 4 - 4 XG 


4X 


2 ZXG 


5X 


(4.7) 


By choosing appropriate combinations of the coefficients cj 1 ' 1 and in the absence of the poten¬ 
tial, one can check that there exist exact linear solutions 0o oc t to (4.3) with F( A"o, Zq) = 0, 
sourcing an exact de Sitter space. Upon turning on a small enough potential for 0, these 
backgrounds can slightly deviate from de Sitter, F[ H 2 . Moreover, if both X Q and Zq are 
of order one on these solutions, the background curvature is of the order H 2 ~ A^/M^ and 


all terms involving (j) in (4.1) contribute equally to the energy density. 


More generally, once shift-symmetry is broken, the parameters X and Z measure the 
departure of quasi de Sitter solutions from slow-roll inflation. In particular, for A" 1 and 
Z 1 one reproduces the standard potential-dominated slow-roll regime. In the opposite 
limit and for small enough V (</>), the acceleration can be supported by the kinetic part of the 


action. At the same time, we have argued below Eq. (3.12) that the same two parameters 


control the magnitude of quantum corrections. Indeed, in the X ^ 1 and Z <?1 regime loop 
corrections are fully under control, which can be seen by estimating typical loop-generated 
terms, e.g. of the form 


(V 2 0) n 


Z n A 


{d<t>) 


2 n 


3 5 


MpjA ^ -4 


X n A 4 


(4.8) 


For X and Z less than or of order one, these are parametrically suppressed with respect 
to the background energy density contributed by the classical Lagrangian, which, as follows 


12 


Note that both G 4X /X and G$x/X in this equation are finite even in the X —> 0 limit due to our 


assumption that G 4 and G 5 start at least quadratic in X , see the discussion below Eq. (3.9). 
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from Eq. (4.6), generically satisfies 3 M^H 2 > XA Saturating either one, or both of the 
bounds in Eq. (3.12) leads to a non-linear regime where various next-to-leading order terms 
in the derivative expansion become important, while quantum corrections are still under 
control. This case will be of our prime interest in the rest of the paper. In particular, in the 
next section we will explore the properties of the corresponding backgrounds in the context 
of inflationary physics. 


5 Inflation 


5.1 The effective field theory for perturbations 


The consequences of approximate galileon symmetry for the physics of the early universe 
can be conveniently studied in the effective held theory of inflation (EFTI) framework [1] [2]. 
Operationally, the EFT of interest is conveniently constructed in the unitary gauge, i.e. by 
choosing the special coordinate system where the time slicing coincides with the uniform 
inhaton hypersurfaces. The spacetime metric can be decomposed into the ADM variables, 


ds 2 = —N 2 dt 2 + 7 ij ( N l dt + dx 1 ) (N^dt + dad) , 


(5.1) 


where N is the laspe and N l the shift, while 7 ^ is the 3d metric. 

The ‘breaking’ of time diffs allows to include objects that are not invariant under 
the full 4d diffs, but preserve the reparametrizations of the 3d coordinates on equal-time 
hypersurfaces. Thus, the possible building blocks of the EFT of interest are homogeneous 
functions of time, intrinsic and extrinsic curvatures of equal-time hypersurfaces, respectively 
^Rij and K.jj , and temporal components of tensors with upper indices such as N = 1 /\J — f/ )() 
na. To construct the action up to quadratic order in metric fluctuations, we perturb these 
quantities around a flat FRW background, ds 2 = — dt 2 + a 2 (t)dx 2 } by writing 5N = N — 1, 
5Kij = Kij ~ H'lij (K = K\), while ^R tj is already at least linear in perturbations. In 
terms of these quantities, the action reads mmmm 


S= / d 4 x^N 


M 2 1 


/M ( (3) fl + K ,3 K,y - K 2 ) - 2/(i)4 + ^ - A(f) 


+ 4Lhbjv 2 - M 3 (t)SK5N - ^-A-(SK 2 - SK^SKtj) + ( 3 )rsn 


2 

M r2 {t ) 


(5.2) 


(6K‘ 2 + SK lJ SK tj ) + m x {t)^R5K + ... 


The ellipses stand for Lagrangian terms constructed of space or time derivatives of N , ^Rij 
and K^. Moreover, we have omitted the operator ^Rij5K l \ since at the quadratic order it 
can be rewritten as a combination of other operators in (5.2), see 
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As far as inflation is concerned, one can always remove the time-dependent function 
/ by means of a conformal transformation of the metric [2j. Thus, in this section we will 
assume / = 1, in which case the time-dependent coefficients c and V in the first line - those 
that lead to tadpoles - are uniquely determined in terms of the background equations of 
motion, 

c = -M^H , A = M^(3H 2 + H) . (5.3) 

The rest of the action contains operators whose coefficients are in principle unconstrained. 
Depending on which subset of these coefficients is nonzero, the above effective theory de¬ 
scribes the linear perturbations in different (and all possible) models of single-field inflation. 

Given that a generic model of inflation can be defined via the set of EFT coefficients 
M 4 (t ), M 3 (t) , M 2 (t) ,..., it is important to understand what are, if any, the possible hier¬ 
archies among these. Indeed, the basic properties of the inflationary background and per¬ 
turbations such as stability, power spectra, non-gaussianities, etc., crucially depend on the 
relative magnitudes of these terms. In the case that the problem involves widely separated 
scales (such as, for instance, the Hubble scale and the UV cutoff of the theory), it is not a 
priori clear what the quantum-mechanically stable values of these Wilson coefficients are. 
In particular, treating them as arbitrary parameters can easily lead to results that require 
tuning, once quantum corrections are taken into account. 


The easiest way to arrive at various hierarchies in quantum field theory is using sym¬ 
metries - exact or approximate. For example, chiral symmetry, although inexact, is the 
central reason behind the technical naturalness of fermion masses in the standard model of 
particle physics. In the same spirit, one can use various symmetries in the effective theory 
(5.2) to constrain the EFT coefficients in different models of inflation. One obvious possible 
symmetry are constant time translations, t —> t + c, which would constrain all coefficients 
(including the Hubble rate) in the EFT to be time-independent, enforcing the background 
solution to be an exact de Sitter space. Another symmetry one can impose are arbitrary 
space-independent reparametrizations of time, t —> /(£); see, e.g., [25]. In the exact symme¬ 
try limit, this would require M, M and fh to vanish, leading to interesting theories in one 
corner of the inflationary theory space [26]. On the other hand, even if time reparametriza¬ 
tions are not an exact symmetry, their weak breaking (appropriately defined in what follows) 
would guarantee that any small values of these coefficients are technically natural. 


Building on the insights of the previous sections, we wish to investigate the conse¬ 
quences of yet another possible approximate invariance of the action (5.2) - that induced 
by internal galileon transformations realized on the foliation scalar as in (1.1). If this scalar 
acquires a linear profile </> oc t driving a de Sitter background (similarly to what we discussed 
in the previous section), the galileon symmetry (1.1) will have the following manifestation 
in the unitary gauge 


t ->■ t + b^x 11 , (5.4) 

where b M and c are a constant four-vector and a scalar, respectively. For more complicated <fi- 
prohles, the associated unitary gauge invariance will have a more complicated form. However, 


14 








in practice it always makes life easier to restore the foliation scalar in the effective action 
(5.2) in order to explore its approximate symmetries. We will take this route in what follows. 


5.2 ‘Ordinary’ effective theories 


Before turning to the case of models with galileon symmetry, it is instructive to recall what 
happens in more standard scalar effective theories. There, one starts with an action organized 
in derivative expansion containing a single scale A - the EFT cutoff, 


S 


Seh + / d A xyJ^g 


A {X) - E(0) + 


w 


+ 


(5.5) 


where, as before, G 2 is a dimensionless function of the variable A" = — (<90) 2 /A], parametriz¬ 
ing the leading derivative effects. 


As discussed in the previous section, for a flat enough potential the scalar background 
profile can be approximated by a linear function of time 0 O = c(t)A 2 t, where c(t) is a slowly 
varying function. This is true both for slow-roll inflation, where dhuf)/dN ~ dlnH/dN is of 
the order of the slow-roll parameters, and for models such as the ghost condensate, where the 
held ‘velocity’ is explicitly constant, i.e. c = const, at the leading order [27] , For definiteness, 
let us concentrate on the latter class of models and assume that c ~ 1 solves G 2 x(c 2 ) = 0. 
On such solutions, the Hubble rate can be estimated as Mp X H 2 ~ A] and is completely 
determined by the leading G 2 (X) piece, higher-derivative operators being unimportant for 
the backgrounds at hand. For example, the quadratic operator shown in Eq. (5.5) can be 
estimated as (□</>o) 2 /A 2 ~ A^(H/Mp\) and thus it is negligible for H <C Mp\, as required for 
a consistent classical description. 

Restricting to the unitary gauge, S(p{t, x ) = 0, one straightforwardly finds that all EFT 
coefficients are determined by the cutoff of the theory 


M 4 ~ 


A 4 


M 3 ~ A 3 


M 2 ~ 


A 


M' 2 ~ 


A 


etc. 


(5.6) 


In this case, the dynamics of small perturbations is fully dominated by the only zero- 
derivative quadratic operator in the effective theory - SN 2 - and most of the phenomenol¬ 
ogy is thus determined by the single coefficient M 4 . This leads to interesting characteris¬ 
tic features, such as the possibility of small speed of sound c 2 of scalar perturbations (if 
|M 4 | Mpj H |) and, associated to it, large non-gaussianities /nl ~ 1/c 2 [2H [2j. The 

contribution of higher-order terms amounts to only slightly correcting the leading results. 
For example, one can show that the correction to the speed of sound from the operator M 2 
(which arises entirely due to mixing with gravity) is of order Sc 2 ~ M 2 /Mp l ~ H/Mp\ and 
can be safely ignored. Consequently, one can consistently consider the perturbations of the 
inflaton held as weakly coupled over a sufficiently broad range of distances encompassing the 
Hubble scale and straightforwardly apply the derivative expansion. This would correspond 
to the case of /^-inflation [29], ghost inflation ra and other related models. 
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£WBG 

M 4 A G 

M 3 HAfl 

M 2 H 2 A^ 4 

m 2 H 2 A 2 4 

1 = 2 

4:X' 2 G2xx 

X 

X 

X 

1 — -H 

II 

OO 

-ex Z ( g 3X + 2 XGsxx) 

-2 XZG zx 

X 

X 

I = 4 

24A ~Z 2 ( 3 G 4 XX + 2A Gixxx) 

SAZ 2 + 2 G 4 xx^j 

— 4:Z"G4X 

M 2 H 2 Afl 

1 = 5 

2 A Z ^3 y + 12G.5 aw + 4A G^xxx^j 

2XZ 3 + 2 G 5X x) 

—2Z 3 G$x 

M 2 H 2 Afl 


Table 1: Contribution from the Lagrangian terms £) VBG to various unitary-gauge operators, defined 
in Eq. (5.2). The relation M 2 = rh 2 holds for all of these. The subscript X means differentiation 
with respect to X and all derivatives are evaluated on the background. We have assumed that all 
background quantities obey ‘slow-roll’, i.e. d/dt <C H, so that terms involving derivatives of X and 
Z have been neglected. Moreover, generically contributes a weakly time-dependent piece 

to the effective Planck mass through the function /(t) in Eq. (5.2). According to our assumption 
about the form of G 4 discussed in Sec. [ 3 J this is at least of order X 2 for small X , so that its 
contribution can be neglected. An analogous piece coming from is suppressed by a further 

factor of the slow-roll parameter. 


5.3 Inflation with WBG symmetry 


Going back to the theories characterized by WBG invariance (4.1), one can straightforwardly 


find the values of the various coefficients in Eq. (5.2), by means of reformulating the discussion 

The contribution from each 


of the previous section into the effective theory language 
Lagrangian Cy BG to different EFTI operators is reported in Tab. [lj 

Interestingly, as a generic property of Horndeski theories, only one combination of the 
operators 5K 2 , 8K l3 8K 13 and ( 3 >R5N, given by the following expression 


5K 2 + 5K ij 5Kij + (3) RSN , 


(5.7) 


appears in the EFT action [30]. Thus, the parameters M 2 and rh 2 are related to each other 
as M 2 = ?7i 2 , while all the others are zero, M’ 2 = = . .. = 0. Moreover, it follows from 

the discussion of Secs. [2] and [3] that this ‘tuning’ is perfectly stable under loop corrections, 
as far as A ^ 1 and Z ^ 1. This implies, for instance, that M 2 M’ 2 can hold in the 
full quantum theory due to the (weakly broken) galileon invariance. Moreover, the unique 


combination (5.7) is in fact a redundant operator and can be removed by a redefinition of the 
metric isniEJ0 As a consequence, the single operator 5N5K - associated with the cubic 
galileon - is responsible for all the differences with respect to the more familiar DBI/ghost 
inflation models, as far as stability and power spectra are concerned. 


13 This redefinition corresponds to transforming into the frame where the tensor modes propagate at the 
unit speed of sound [32]. 
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With this in mind, we derive the quadratic action for the comoving curvature perturba¬ 
tion ( (defined in the standard way, 7 ^ = a 2 e 2< ^8ij) for the general case of nonzero coefficients 
M 4 and M 3 in Appendix [bJ The full result (away from any decoupling/short-distance limit) 
reads 


where 


S(= d 4 x a 3 


A(t) C 2 - Bit) 


(VC ) 2 


Mpj(—4Mpji/ + 3 M 6 + 2Mp l M 4 ) 

(M 3 - 2M 2 l H ) 2 

Mpj(—4Mpji/ + 2Mp l HM 3 - M e + 2M 2 l d t M 3 ) 
(M 3 - 2M 2 l H) 2 


(5.8) 


(5.9) 

(5.10) 


and the speed of sound for short-wavelength perturbations is simply c 2 = B/A. Depending 
on the precise values of the parameters A" and Z, there are two phenomenologically distinct 
regimes in a theory of inflation with WBG invariance, which we review in what follows. 


5.3.1 The kinetically driven phase 


The strongly coupled backgrounds - i.e. the ones for which all terms (3.1)—(3.4) are of the 
same order - correspond to the function F{ A"o, Zq) vanishing at the leading order, as already 
noted above. For X and Z strictly constant, one has an exact de Sitter space, which can 
be made quasi-de Sitter by turning on a small potential for 0. One can see from Eq. (4.4) 
that for F to vanish requires that at least Zq be order-one. Although less evident, it is 
easy to show that also A" 0 has to be order-one, if the potential is to provide a sub-leading 
contribution to the energy density. Indeed, it follows from Eq. (4.6) that the energy density 
scales as Mp Y H 2 ~ h.\X in this case. Together with the following equality, 


z = Vx 


Mp\H 

~w~ 


(5.11) 


which simply follows from the very definitions of A" and Z, this implies that Z ~ A". One 
thus concludes that whenever the theory approaches strong coupling on a (quasi) de Sitter 
background, Z ~ 1 , all functions Gi should be generically resummed. We stress again that 
this only concerns the case when V ( 0 ) can be neglected at the leading order - a condition 
we will give up below. 


For large backgrounds that nearly saturate the bound (3.12), one consequence of the 
approximate symmetry under ( 1 . 1 ) is the following radiatively stable hierarchy 


M^H 2 ~ M 4 ~ M 3 H ~ M 2 H 2 ~ m 2 H 2 » mi H 3 , M' 2 H 2 , ... , (5.12) 
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as can be seen from Tab. [TJ Here, ellipses denote all other coefficients of operators of 
various canonical dimension, properly normalized with the help of the appropriate number 
of powers of the Hubble scale. The strongly-coupled case at hand corresponds to the largest 
possible values of the coefficients M 4 , M 3 and {M 2 ,m 2 }, i.e. of order Ad^H 2 , M^H and 
Mpj respectively. For the latter values of these operators, their effects become order-one 
important for physical quantities such as, for instance, the power spectrum or the speed 
of sound of the scalar perturbations - see the exact expressions (5.9) and (5.10). In this 


case, derivative interactions either dominate or contribute significantly to the background 
curvature, unlike potentially-driven models such as for instance slow-roll or DBI inflation. 
Remarkably enough, for a sufficiently large hierarchy between A 2 and A 3 - or, equivalently, 
between M pi and A 3 , the strong coupling of our backgrounds does not necessarily imply 


the breakdown of the classical description, as follows from the estimates of Eq. (4.8). For 


example, in terms of the EFT language, the operator (□</>) 2 /A 2 leads to the 5K 2 term of 
order 5M 2 At/A! Mpj(if/Mpi) 2//3 , which is suppressed by a tiny factor compared to the 
leading contribut ioi| T ~ i ] 

There is another important characteristic to the ‘large’ backgrounds X ~ 1, Z ~ 1 


- or, equivalently, to large values of the EFT coefficients given in Eq. (5.12). Assuming 
MpjFT ^ M 4 , as required by being close to de Sitter space, there is no short-distance limit in 
which the metric can be decoupled from the inflaton perturbations. Therefore, mixing with 
gravity is order-one important at all scales. To see this, let us restore the Goldstone boson 
7r, nonlincarly realizing the spontaneously broken time diffs: 


g 00 ->■ g 00 = -N~' z = -1 - 2tt - 7 r z + 


7—2 


•2 


(dn) 2 


d 2 7T 

5K ^5K -- + 0( it 2 ) . 


The relevant part of the n action then reads 


+ \m\6NY - m 3 snsk -> *1 + + • •. . 


(5.13) 


where we have defined the ‘decay constant’ f 2 = (M 4 /2 — M^H) 1 ^ 2 and canonically nor¬ 
malized the Goldstone and the lapse variable as 7r = 7r c // 2 and 5N = 5N C /Mpi, respectively. 
One can see from the last term in this expression, that the mixing with gravity is indeed 
important for the Goldstone dynamics at all distance scales if the EFT coefficients saturate 
the strong coupling bound (i.e. if M 3 ~ MpjM 2 ). This means that, in this case, one has to 
perform the full analysis including dynamical gravity in order to extract reliable predictions 
from the theory. That gravity is order-one important for such backgrounds can also be seen 
from the exact expressions for the kinetic coefficients of the curvature perturbation (5.9) and 


14 Note that quantum corrections generically introduce higher time derivatives (e.g. SN 2 ) at the same 
scales as the higher space derivatives. These lead to loss of unitarity associated with ghost excitations with 
masses above the EFT cutoff (cured by whatever completes the theory in the UV). The scales associated 
with the higher space- and time-derivative operators can only be disentangled when the speed of sound is 
small, c s < 1, in which case the latter scale becomes enhanced by powers of cj 1 . 
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(5.10). The calculation of the cubic action/non-gaussianities in this regime is much more 
involved due to the non-decoupling of gravitational effects and will be reported on elsewhere. 


In principle, Z may end up being somewhat smaller than one due to an accident/fine 
tuning (as discussed above, in our case any tuning of coefficients is respected by quantum 
corrections in the theories at hand). In this case, one can then define the decoupling limit 
governed solely by the dynamics of the Goldstone mode or, equivalently, by the comoving 
curvature perturbation ( = —Htt. In this case the analysis beyond the quadratic action 
becomes significantly simpler. The non-gaussianity in the decoupling-limit has been studied 
in the context of galileon inflation in Ref. [13] and we will not reproduce the results here. 
Unlike DBI/ghost inflation, this limit contains cubic operators with four derivatives, such 
as for instance <9 2 7r(<97r) 2 . However, all of these operators are redundant and can be reduced 
on-shell to those with three derivatives [33], i.e. n^dn) 2 and n 3 . This means that, at least 
in the decoupling limit, galileon inflation does not lead to new shapes of non-gaussianity, 
distinct from those that already arise in DBI models. The genuine novelty is that, unlike 
DBI theories, more than one operator becomes nonlinear in galileon inflation, leading to a 
richer parameter space. As a result, the amplitude of equilateral non-gaussianity //jj[ ul is not 
necessarily linked to the speed of sound of scalar perturbations and it can be larger than the 
standard DBI result, /l// 111 rs_/ 1 /c 2 , see [I3| for a more detailed discussion. 


5.3.2 Slow-roll inflation with WBG symmetry 


A different and particularly interesting case corresponds to a potentially driven universe with 
a subleading contribution to the energy density due to the derivative part of the action, 

Ml x H 2 ~ V > . (5.14) 


In the inflationary context this is just the slow-roll limit. The slow-roll parameter e can be 


estimated from the Friedmann equations (4.6) and (4.7) to be of order 


H 




A‘A' 


//■ 


(5.15) 


where we have made use of the fact that F is an order-one function and have assumed that 
F and </> vary slowly with time, F <C HF, f <C Hf. Using the identity (5.11), the last 


expression for e immediately implies that Z has to be parametrically larger than X, 

Ve 


(5.16) 


Let us consider now a potentially-driven regime, characterized by Z ~ 1. According to 


Eq. (5.16), A" is parametrically suppressed, being of order -y/e, and it follows from (5.15), 
that 


Mp } H\ = Mp Y H 2 e ~ A 4 2 X . 


(5.17) 
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The latter scale, denoted by f n in Eq. (5.13), defines the scalar two-point function in ordinary 
slow-roll models, setting the normalization for the Goldstone mode of broken time diffs, 


C = + ... , 


(5.18) 


where ellipses denote all sorts of additional terms, completely unimportant in slow-roll in¬ 
flation, including those describing the mixing of tt with gravity. 


Let us now have a look at the values of the additional operators present in the case 
of interest and listed in Tab. [l] A crucial observation is that, when Z is of order one, the 
contribution to M 4 and M 3 from all jCY BG with 7^2 scales as M 4 ~ A 4 A" ~ M^H 
and M 3 ~ Mp X H/H respectively. The M 4 operator that arises from £^ BG is down by 
an additional factor of X ~ i/e compared to this, so that, as long as perturbations are 
concerned, higher-derivative terms are more important than the next-to-leading order terms 
in the derivative expansion of the form (<90) 2n . Such a scaling of extra EFT operators is quite 
remarkable, since it implies that the perturbation Lagrangian (5.18) is modified at order one 
compared to the usual slow-roll case, while the background evolution is still fully governed 
by the potential 3 ) 


This can also be seen from the exact expressions for the quadratic scalar perturbation 
Lagrangian. Taking into account the magnitudes of the two EFTI operators at hand, let us 
denote 


M 4 = —aMp\H = aM^H 2 e , 


TT 

M 3 = 1 3M 2j— = -/3M 2 } He , 

H 


(5.19) 


where a and (3 are constants of order one. Then, from Eqs. (5.9) and (5.10) one finds 

(1 + !)<«?,+ e>( £ 2 ), s=( i-0 


A = 


- eM 2 l + 0(e 2 ), 


r 2 = lzl 

s 2 + ex 


(5.20) 


This explicitly shows that, by suitably choosing the order-one coefficients a and /3, which 
themselves are functions of the free parameters of the theory, the speed of sound of the scalar 
perturbations end up anywhere in the allowed region 0 < c 2 < 1. Moreover, tuning f3 to be 
close to 2 yields a regime with highly subluminal scalar perturbations, potentially implying 
strong coupling and large ‘accidental’ non-gaussianities pH], all in a technically natural way 
(recall that any tuning in the theories at hand is protected by WBG invariance). 

Last but not least, we note that in the potentially driven regime, one can generically 
define the decoupling limit - unlike the kinetically dominated case considered above. This 
can be seen by noting that the term describing the mixing of the (canonically normalized) 
Goldstone tt c with gravity is parametrically suppressed on backgrounds close to de Sitter. 

15 This is analogous to what occurs in DBI inflation m , even though there are crucial differences. Here 
we have in mind a scenario closer to slow-roll inflation than DBI models, in that e.g. the Hubble friction is 
important for the inflationary dynamics of <fi and sustaining inflation requires a flat enough potential - or, 
in other words, that the usual relation e ~ (Mp\V'/V) 2 holds. This is not necessarily true in DBI inflation. 
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Indeed, using the expressions for the magnitudes of EFT coefficients in Eq. (5.19), one can 


readily evaluate the coefficient of the mixing term in Eq. (5.13) to be of order 


M 3 


M^fl 


V'i 


(5.21) 


During inflation, this is a parametrically small number meaning that the Goldstone mode can 
be studied independently from the gravitational degrees of freedom, at frequencies of order 
of the Hubble scale. The same conclusion can be made regarding all other EFT operators 
present in the theory. 


One can also see from Table [T| that, while generalizing galileon inflation, the most 
general inflationary theory defined by ( |4.1| ) possesses a much broader parameter space. While 
three parameters in Tab. [I] characterize galileon inflation {Gzx, G 4 xx, G$xx), a generic 
theory with WBG invariance has six more parameter^*] at the level of the two-point scalar 
function (G 4 x, G 5 x, G2xx, G 3 xx, G 4 xxx, G 5 xxx)- Higher-point functions will introduce 
even more freedom. Moreover, arbitrary tuning of the parameters is stable under radiative 
corrections and taking advantage of this fact can lead to interesting theories in various corners 
of the EFT of inflation theory space. The complete analysis of possible physical predictions 
of the theory (4.1) is beyond the scope of the present discussion. 


6 Late-time universe 


6.1 Dark energy 


The theories characterized by WBG symmetry can also be applied in the context of the 
late-time cosmic acceleration. For simplicity, we will neglect the scalar potential V (0) in 
this section, contrarily to the case discussed in the previous section. Moreover, we reabsorb 
the canonical scalar kinetic term in the definition of the function Gb, so that the action of 
interest reads as follows 


S 



1 

2 


1=2 


( 6 . 1 ) 


The homogeneous equations on a flat FRW background obtained by varying this action are 
similar to those, discussed in Sec. [5] However, in this section we choose to rewrite them in 
a slightly different form and notation. 


In terms of the two variables X and Z defined in Eqs. (3.5) and (4.5), the homogeneous 
scalar evolution equation reads 

d 


dt L 


^0F DE (X, Z) 


= 0, 


( 6 . 2 ) 
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Assuming G 4 x and G are zero, which we did in this paper. 
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where i+> E (X, Y) is defined as 

+ ZGsxx'j • (6.3) 

We assume that all matter couples minimally to the metric g iwi so that the Friedmann 
equations can be written as 


F de (X, Z) = G 2X - 3 ZG iX + QZ 2 


G 


4X 


X 


+ 2 G 


4XX 


+ Z d 


G 


5X 


X 


3 M 2 H 2 — p m + pde , (6.4) 

-2 M 2 H = Pm + Pm + p DE + p DE , (6.5) 


where p m and p m respectively denote the energy density and pressure of matter. The effective 
Planck mass squared, M 2 , is generically time-dependent due to the non-minimal couplings 
of the scalar in £^ BG and £^ BG , and is explicitly given by the following expression 

M°; = M El (1 + 2 G 4 - 4 XG 4X - 2 ZXG 5X ) . (6.6) 


In terms of the EFT parameters of Eq. (5.2), Ml = M^f + M 2 . Finally, the unperturbed 
effective energy density and pressure read 


Pde — 


+12 Z' 1 


X + 2 a 2X 


6 ZG 


3X 


G. 


4X 


X 


+ 2 G 


4XX 


+ AZ d 


G 


5X 


X 


+ G 


5XX 


Pde + Pde — 2 A 2 XFDE + 2M E \(j)X 


x 


G 3 .y - 4Z 


G 


4X 


X 


+ 2 G. 


4XX 


-Z 2 


.G 


5X 


X 


+ 2 G 


5XX 


(6.7) 


( 6 . 8 ) 


For Horndeski theories, the phenomenological deviations from the ACDM model at 
the level of linear perturbations can be conveniently parametrized in terms of four time- 
dependent parameters [22] • It is convenient to use the dimensionless functions introduced in 


, which can be expressed in terms of the EFT coefficients of Eq. (5.2) as 


a.K = 


2c + M 4 


OiB — 


M 2 J - M 3 


OiM = 


M 2 J + M 2 


M 2 


CtT — 


(6.9) 


M|ff 2 ’ - 2 m;h m?h a/; 

where the function c is fixed by the background Friedmann equations derived from action 


(5.2), 


2c — pde + Pde + M El (Hf — f) . (6.10) 

For the theory of interest (6.1), their explicit expressions are reported in App. [C] The con¬ 
tributions from the Lagrangians jCY bg to each of these parameters is given in Tab. [ 2 ] 

Let us first briefly comment on the physical meaning of the ^-parameters and their 
observational consequences (see, e.g. [36, 30] for a more detailed discussion). 
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£WBG 

M t 2 H 2 A^ 4 a K 

M 2 H 2 A 2 - 4 a fl 

M 2 H 2 A^a M 

M 2 H 2 A 2 " 4 a T 

1 = 2 

2 X(G 2X + 2 XG 2XX ) 

X 

X 

X 

I = 3 

-12 XZ(G 3X + XGzxx) 

XZG 3X 

X 

X 

I = 4 

12 XZ 2 ^ — + 8G4XX + / ^XG\xxx'\ 

—4XZ 2 + 2G 4X x^j 

8 

M B 

4Z 2 G ix 

1 = 5 

4XZ 3 (3^ + 7 G 5XX + 2 XGsxxx) 

-XZ 3 (3^ + 2 G S xx) 

9 HZ* i 

2 H 2 G 5 X 2 A 3 a B 

2 Z 2 ( z ~m) Gsx 


Table 2: Contribution from the Lagrangian terms £) VBG to the various effective a-parameters 
defined in Eq. (6.9), denoting the deviations from the ACDM model in linear perturbation theory. 
Their complete expressions are reported in App. |C| see Eqs. (C. 1)—(C.4). 


• The first of these functions, aix, parametrizes the kinetic energy of scalar perturbations, 
induced by the four Lagrangian terms CY BG ■ I 11 terms of the EFT for perturbations 

), the contributions to olk arise from the operators oc 8 N 2 , described by the co¬ 
efficients c and M 4 . This parameter is enough to describe linear perturbations in 
the minimally coupled quintessence and ^-essence models. In the minimal case, dark 
energy fluctuations behave as those of a perfect fluid with the speed of propagation 
c 2 = 3(1 + Wde^de/ciag where IIde is the energy density of (j) in the units of the 
critical one. A well-studied example corresponds to the limit olk 3> 1, which leads to 
zero sound speed PEZ]. 

• The second function, a B , parametrizes the mixing between metric and scalar field 

fluctuations [38] . It can be induced by the Lagrangians with 3 < / < 5, or by 

the operator containing 8 N 8 K in the EFT for perturbations, which describes kinetic 
mixing with gravity. This operator induces a typical scale (approximately given by 
k B /a ~ 3i/^/f2 m /2 for % = = 0) [35], below which dark energy can cluster with 

energy density fluctuations dpoE — 2 8p m a 2 B / ({ol 2 k + 6a B )c 2 ). Above this scale, dark 
energy fluctuations simply behave as those of a perfect fluid. 

• The function is defined as a m = d In M 2 j In a, and is thus related to the time 
variation of the effective Planck mass M 2 induced by the Lagrangians £^ BG and 

It parametrizes the non-minimal coupling in scalar-tensor theories such as Brans-Dicke 
[39] and, as such, induces a slip between the gravitational potentials. In terms of the 
EFT operators, olm can be generated by time variation of the parameters / and M 2 . 
In particular, in f(R) theories, ckm = — ol b . 

• Finally, the last function, parametrizes the deviation of the speed of sound of tensor 
perturbations from unity in the presence of the Lagrangian BG and £^ BG [35] , and 
also induces a slip between the gravitational potentials. It is generated by the EFT 
operator, multiplied by the coefficient M 2 . 
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A particularly interesting case corresponds to a constant X, where the scalar held 
grows linearly with time, (fro oc t. In this case, the functions Gi, and all their derivatives are 
time-independent when evaluated on the background, while Z changes proportionally to the 
Hubble rate. One can see from the form of the scalar equation, Eq. (6.2), that the 0o oc t 
profile can only be a solution if the following relations are satisfied 


G 2 x — 0 


G%x — 0 


Gax + 2A"G' 4 va' — 0 


3 G 5 x + 2 XG 5 xx — 0 . (6.11) 


It then follows from Eqs. ( |6.7 ) and ( 6 . 8 ) that the effective energy density on such backgrounds 
is given by /?de = —A^(G 2 + ^Z z G^x) and the equation of state is that of the cosmological 
constant, with wde = Pde/pde = —1- Moreover, while the dimensionless parameters «a', 
a.M and ax in Tab. [2] are generally nonzero and time dependent, Eq. (6.11) implies that 
as = 0. In this case, in order to avoid ghost instabilities one must require that ax > 0 
[231136]. Note that since Z oc H, pde is time dependent, even though wide = —1- This is 
because the effective Planck mass varies with time and the dark energy does not follow the 
usual conservation equation. Demanding that the energy density remains finite at early times 
requires it to be constant, i.e. G 5 x = 0 and thus that also = 0, see Tab. [2} I 11 this case, at 
the background level this model exaclty behaves as a cosmological constant. At the level of 
linear perturbations, the dimensionless parameters ax and ax do not vanish and contribute 
to the sound speed of scalar and tensor fluctuations, respectively as c 2 S = — 2 ax/ax and 
Cx = 1 + «t (see [36| [30]). One can avoid gradient instabilities by requiring the positivity of 
both these speed of fluctuations squared, which implies — 1 < ax < 0. 

This is the simplest application of WBG symmetry to the late time-acceleration. The 
assumption that the background profile of the held is linear with time, 0o cx t, and that the 
background expansion history is the same as in ACDM considerably restricts the values that 
the parameters a a can consistently take. These solutions are not possible for the covariant 
galileon [2]. In this case one usually assumes a tracker solution with Z = const. [TO]; 
this imposes a particular expansion history, which has been shown to lead to observations 
that are disfavoured with respect to ACDM (see m for a recent analysis and e.g. @2103] 
for analysis that do not assume the tracker solution). More sophisticated examples can be 


constructed using the Lagrangians (3.1)—(3.4), by allowing the background solution for (j) to 


be different from a profile linear in time. In this case, from Tab. [2] one generically expects 
that the parameters a a may assume any value smaller than unity, a a ~ X 1. 


6.2 Static sources and the Vainshtein mechanism 

Any theory with gravitationally coupled scalars should be able to pass the solar system 
constraints in order to be able to qualify as an acceptable alternative to general relativity. 
In the context of theories with WBG invariance, one can imagine coupling (j) to the trace of 
the matter stress tensor in the following way 

Grnat ~ wj—0T m . (6.12) 

M Pi 
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It is then necessary to make sure that the exchange of 0 by realistic sources does not modify 
the Newtonian potential by any detectable amount. There are a number of ways to suppress 
the contribution from 0 to the gravitational potential of astrophysical objects (stars, planets, 
etc.) in modified gravity. In the case of the theories with approximate galileon invariance 
we are interested in, what guarantees that 0 is screened beyond the observable values is the 
Vainshtein mechanism [33], 35 ]- Indeed, in the decoupling limit (3.8), the theory reduces to 


the Galileon, which is known to naturally incorporate the screening below the Vainshtein 


redius pj. For a souce of mass M S( 


this radius reads 


rv = 


M s< 


Me, J 


\ 1/3 

- I 


(6.13) 


The decoupling limit, on the other hand, does capture all relevant astrophysical scales and 
all results obtained in that limit can be fully trusted. To see this directly, one can evaluate 
the quantity X on a pure Galileon solution of Ref. Pj and check that it is extremely small 
everywhere in spac^J so that all the results of Ref. pj apply to our case without any 
modification. 


7 Conclusions and outlook 

Galileon-symmetric scalar theories have played an important role in cosmological model¬ 
building in the context of both the early and the late-time universe. The exact symmetry 
is only possible in the absence of gravity and it leads to remarkable properties - both at 
the classical and quantum levels. However, upon inclusion of gravity galileon invariance 
necessarily has to be broken, raising the question regarding the fate of its attractive flat- 
space properties. 

In this paper, we have introduced the notion of weakly broken galileon invariance , which 
allows to naturally generalize the galileon to curved space (see Sec. [2]). Requiring that the 
quantum non-renormalization properties of the flat-space theory is preserved to the maximal 
possible extent also in the presence of gravity has led us, quite uniquely, to a particular sub¬ 
class of the ghost-free Horndeski theories, presented in Sec. [3] The resulting effective theory 
is characterized by two scales - or, more appropriately, by parametrically small coefficients 
of symmetry-breaking operators in the units of the effective held theory cutoff - allowing 
to retain the quantum non-renormalization properties of the galileon for a broad range of 
physical backgrounds. 

As we have shown, the latter properties guarantee that de Sitter solutions that these 
theories generically possess are insensitive to loop corrections. This gives a lot of freedom to 

17 In the case of the cubic Galileon, X becomes of order one only at the Schwarzschild radius of the source 
rs, while, in the most general case, it is a constant within the Vainshtein radius, of order of the tiny ratio 
(rs-ffo) 2 / 3 , where H 0 is the current Hubble rate. 
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implement these solutions for building technically natural models of both the inflationary, 
as well as the late-time universe. One of the purposes of this work has been to introduce 
these models and set the stage for their detailed phenomenological studies. 

In Sec. [5j we have studied two relevant regimes in the context of inflation. The first 
one describes genuinely strongly coupled physics, characterized by A" = </> 2 /A| ~ 1 and 
Z = Hcf)/ A| ~ 1 in the notation of Sec. [ 3 } It generalizes galileon inflation of Ref. [T3] and 
saturates the bound on the ‘strength’ of the backgrounds for which the classical analysis 
can still be trusted. The dynamics of perturbations on these backgrounds is completely 
dominated by higher-derivative operators. Yet, the mechanism is fully predictive, quantum 
loops giving subleading corrections to the classical results. I 11 a generic theory with weakly 
broken galileon invariance, there are six additional parameters (two of which have been set 
to zero in this paper) compared to the particular case of galileon inflation at the level of the 
scalar two-point function. More freedom arises when higher n-point functions are concerned. 


The other, particularly interesting, regime corresponds to partially strongly coupled 
backgrounds, with Z ~ 1 and X ~ i/e, resembling slow-roll inflation. In this scenario, 
inflation happens mostly due to the potential; however, weakly broken galileon invariance 
makes it possible for the higher-derivative operators to become relevant for the dynamics 
of scalar perturbations, modifying their action at order one. In particular, this allows for 
a significantly subluminal speed of sound for the curvature perturbation ( as a result of 
a possible parameter tuning - which is, importantly, respected by loop corrections to the 
leading order in Mfl 1 . The small speed of sound generically translates into strong coupling of 
scalar perturbations, leading to large non-gaussianities. Such ‘accidentally’ moderate/strong 
non-gaussian perturbations can never arise in ordinary slow roll models. Thus, this represents 
an interesting possibility that the theories with weakly broken galileon invariance can offer. 
These matters will be studied in detail in an upcoming work j3Tj. 

The theories introduced in Sec. 3 can also be used to drive the current accelerated 
expansion of the universe, as discussed in Sec. [6j The situation is richer than the covariant 
galileon case studied in the literature. In particular, while one can recover a background 
expansion history closed to the one currently observed, at the level of linear perturbations 
these theories can lead to modifications of gravity of different types. Expressed in terms 
of the dimensionless effective field theory parameters a a defined for Horndeski theories, 
one generically finds that these cannot be too large, a a ^ 1, leaving the possibility of 
modifications which are compatible with current constraints but possibly testable by future 
observations. We leave a more thorough study of concrete scenarios for the future. 


As a final remark, we note that we have not carried out a complete analysis of quantum 
corrections in the theories with weakly broken galileon invariance; we have rather made the 
most conservative estimates of their magnitudes, so that the statements made above should 
remain true for the most general completion of these theories above the scale A 3 . Whether 
or not our findings regarding the structure of the low-energy effective theory (4.1) imply 
anything useful about its UV completion, is still to be seen. We postpone the detailed study 
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of this question for future work. 
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A Quantum corrections 

In this appendix, we extend the discussion of Sec. [2] to show that the vertices with three 
solid lines and one graviton, as well as five solid lines and two gravitons (see Fig. [ 2 ]) can be 
removed in a suitable curved-space extension of the galileon terms. 

We have considered the case of the cubic galileon in the main text. One has to work a 
little more to understand the case of quartic minimally coupled galileon. The contribution 
to the second vertex from the lower line of Fig. [2] can be extracted by picking up a factor of 
dh from the covariant derivative acting on one of the scalars. The relevant term reads, 


V M V, 


2 (Sfj h i>i/ T d v h ptl dph[iv) d p (f) 


(A.l) 


Here (and in the rest of the present section), by we mean ‘equals up to a total derivative 
and up to terms with fewer factors of d<j>\ and we do not distinguish between upper and lower 


indices for simplicity (everything is contracted with the flat-space metric). Using (A.l), we 
find after a little bit of algebra 


-4£2 


-A(d(f)) 2 h pu dp(f) {d n d v - O/^) <9 p </> , 


(A.2) 


where denotes the Lagrangian term, obtained by minimally covariantizing (2.3). In 
serting the latter vertex into a generic 1PI loop diagram would induce galileon svrnrnetrv- 
breaking operators at the scale Mp{ 6 A^ 6 , parametrically exceeding A 2 - as discussed above 18 
Fortunately, it turns out possible to raise it by adding non-minimal couplings, capable o 
eliminating the vertices with three solid lines coming from the minimally coupled theory. For 
the quartic galileon, the right coupling is y/^g(d<j)) A R . Indeed, expanding the Ricci scalar 
to the linear order in h and integrating by parts, we obtain 


g(dcj)) 4 R ~ 4 (d4>) 2 h pu d p (t) (d p d u - <9 P 0 + Q(h 2 


(A.3) 


18 Although parametrically higher than A 3 , this symmetry-breaking scale is still low enough to disrupt the 


de Sitter spectrum (5.12). 
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implying that the two contributions to the possible vertex with one graviton and three solid 
lines exactly cancel for our resulting generalized theory. 


The way we have chosen to couple the galileon to gravity has been dictated by our 
desire to raise the scale suppressing loop-generated symmetry-breaking operators. However, 
there is something more to it: the resulting non-minimal theory can be recognized as that of 
Horndeski type, leading to second-order equations both for the metric and the scalar. In the 
retrospect, this is not that surprising: the ‘three solid line’ vertices, if present, are the only 
ones that would introduce higher-order equations for the metric in the theory of minimally 


coupled galileon, as can be easily seen from e.g. Eq. (A.2). Our non-minimal term exactly 


cancels this contribution. We thus observe something similar to what happens in massive 
gravity there, eliminating terms that lead to higher-order equations (or, equivalently, 
to the Boulware-Deser ghost) automatically raises the quantum cutoff of the theory; in our 
case, the cutoff (A 3 ) remains unaltered, but a similar procedure raises the scale of breaking 
for the galileon symmetry. 

The same result can be straightforwardly extended to the case of the quintic galileon. 


A brute-force expansion of the minimally-coupled term, using (A. 1 ), yields 


~ (<90) 2 <9 p 0[<9 p h ((D0 ) 2 - (d p d u (j)) 2 ) - 2d p h pu {d p d y (j£2(t) - d p d a (t)d v d a (j)) 
- 2 d p h pp ((D0) 2 - ( d a dpq i>) 2 ) + 4 d^h^^d^d^Dcj) - d p d a (pd v d c 


(A.4) 


It is impossible to manipulate the last expression by partial integration or otherwise, so as to 
be left with less than three factors of <90. This would lead to the last vertex with three solid 
lines in Fig. [2j and therefore higher-order equations for the metric. One can however show 
that this term can be eliminated by adding a non-minimal piece with an exact coefficient 
corresponding to the Horndeski theory. To this end, we note that to the linear order in the 
metric perturbation, the following relation holds 

v/=^(<90) 4 G pi , W0 = -^(00) 4 e M « pA e^ A a Q ^V^^0 + 0(h 2 ) , (A.5) 


where e denotes the totally antisymmetric symbol with 61234 = 1. Expanding the antisym¬ 
metric product and integrating by parts, we obtain 


(00) 4 G^V^V"0 ~ (<90) 2 <9 p 0[<9 p h ((a p ^0 ) 2 - (D0) 2 ) + 2 C0V (<9 P <9„0D0 - c0<9 Q 0c0<9 a 0) ]. 


This exactly cancels the first line of the minimal term’s contribution (A.4). Partially inte¬ 
grating the second line of ( |A.4 ) on the other hand, gives 

- 2(<90) 2 <9 p 0<9„/r pp \d p d a (t)d v d a (t) - d^d^Dcj) + <9 P D0^0 - d^d^d^] . (A. 6 ) 


One can now check, that whatever stands to the right of (<90) 2 <9 p 0 in the last expression, 
is a total derivative (d u of a local operator). This means that one final partial integration 
can get rid of one solid line in the corresponding vertex, reducing the number of solid lines 








to two. This proves our statement regarding elimination of vertices with one graviton and 
more than two solid lines. 


One last statement we wish to prove regards the absence of vertices with two gravitons 
and more than four solid lines. The only possible obstruction comes from the non-minimal 
coupling obtained above, (c^^G^V^V^, in which we take the Einstein tensor linear in the 
metric perturbation, and pick up another factor of h from the covariant derivativ^J One 
can show however, that the corresponding term can always be put into the form in which 
there are no more than four scalars with a single derivative acting on them. Indeed, using 
the expression for Gp„ in terms of the antisymmetric symbol, we have 

~ ^{d^Ye^apxevpaxdadphp^d^h^ - <9 7 /v)<9 7 0 . (A.7) 

This form makes it straightforward to convince oneself, that the factors involving two gravi¬ 
tons in the above expression collect into either a total derivative, or a total derivative up to 
corrections involving at most four factors of d(j)\ 




pX^vPcrxdadfjhppdphryu dp GpopX^v 3 cr A dfj hp a h v ) 


p\^v (3<j\d0.9ph po-dryh^i, ~ -(< 90 ) 4 < 9 7 0 d 1 ( 

£p,a px^-u pcr\dadph 

pcjhfxu) 


(A.8) 


At this point, all that remains is to do a single partial integration to complete the proof of 
the statement, made in the beginning of the present section. 


B Spectrum of (near-) de Sitter space 

Exact de Sitter space saturates the Null Energy Condition (NEC), which seems to be one 
of the most robust ways of covariantly expressing the requirement of energy positivity in 
general relativity]^] This seems to be supported by the fact, that at least in the simplest 
effective theories — those that obey the standard derivative expansion with non-degenerate 
leading-order effects — NEC-violation inevitably implies either ghost or gradient instability 
at short distances |4Uj. Moreover, the exact dS limit of the corresponding backgrounds 
can only be stable over a finite time, necessarily suffering from Jeans-like instabilities [T]. 
Furthermore, there is additional, purely quantum evidence against the existence of exact de 
Sitter spacetime B3- In view of these facts, it is important to understand the most general 
conditions for stability of dS as well as its small deformations. 

In this appendix, we study this question for the effective theory governed by WBG 
invariance. Since these theories lead to second-order equations of motion, the dispersion 

19 A vertex with more than four solid lines coming from the minimally coupled quintic galileon would 
involve at least three gravitons. 

20 This is unlike e.g. the Strong Energy Condition (SEC), which is known to be violated both in the early 
and in the late universe. 
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relation for the scalar perturbations will have a oj 2 oc k 2 form. Moreover, as noted above, the 
operators M 2 and m 2 are redundant, and can be removed by a disformal transformation on 
the tensor mode; one can therefore concentrate on the only remaining quadratic operators: 


8N 2 and 8N8K in the effective theory (5.2). It is convenient to fix the unbroken spatial cliffs 


so as to put the 3D metric in the following form 


9ij ® (^)[(1 "b 2 A T ij\i Ifii 0 j 


(B.l) 


where ( and 7 capture the physical scalar and tensor perturbations. Modes of different 
hclicity do not mix on the homogeneous and isotropic backgrounds considered here, so that 
we can discard the hclicity-1 part of the shift altogether, setting N l = 8 lJ djf8. I 11 the effective 


theory for perturbations (5.2), both the lapse and the shift can be integrated out from their 


respective equations of motion 
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Mi 


pi 


{3) R + K 2 - K ij Kij + —H - 2(3 H 2 + H) 


+ 2 M a 8N - 2 M 3 (8K - 3 H8N) = 0 , 


(B.2) 


Mp! {K l 3 - 8)K) - 8 jM 3 8N - M' 2 8K l 3 - M 2 8K8) 


= 0 , 


(B.3) 


where the lapse perturbation is denoted by 8N = N — 1 and all EFT coefficients except for 
M 4 and M 3 are set to zero for the reasons outlined above. Solving these equations for 8N 


and to linear order and substituting the solutions back into (5.2), one folds the following 
quadratic action for the curvature perturbation 


,5' c = / d^x a 3 


A(t) C 2 - B(t) 


(VC)' 


(B.4) 


where 


+ 3 M 6 + 2MpjM 4 ) 

(M 3 - 2 M 2 X H) 2 

+ 2 M 2 X HM 3 - M 6 + 2M 2 l d t M 3 ) 
(M 3 - 2Mp l H) 2 


(B.5) 

(B. 6 ) 


The speed of sound of short-wavelength scalar perturbations is given by c 2 = B/A. The last 
two expressions agree with those of Ref. [|9] up to a change of operator basis. 

21 For the purposes of deriving the quadratic action for £ only the linearized version of these equations 
matters. This is because the contribution of higher-order terms in SN and (3 to the free action for the 
curvature perturbation vanishes on-shell. 
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C Effective parameters for dark energy 


Deviations from ACDM for WBG invariant theories can be parametrized in terms of the 
following time-dependent functions [36 Fd 
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